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Abstract. The most recent results from the Boomerang, Maxima, DASI, CBI and VSA CMB 
experiments significantly increase the case for accelerated expansion in the early universe (the 
inflationary paradigm) and at the current epoch (dark energy dominance). This is especially so when 
combined with data on high redshift supernovae (SN1) and large scale structure (LSS), encoding 
information from local cluster abundances, galaxy clustering, and gravitational lensing. There are 
"7 pillars of Inflation" that can be shown with the CMB probe, and at least 5, and possibly 6, of 
these have already been demonstrated in the CMB data: (1) the effects of a large scale gravitational 
potential, demonstrated with COBE/DMR in 1992-96; (2) acoustic peaks/dips in the angular power 
spectrum of the radiation, which tell about the geometry of the Universe, with the large first peak 
convincingly shown with Boomerang and Maxima data in 2000, a multiple peak/dip pattern shown 
in data from Boomerang and DASI (2nd, 3rd peaks, first and 2nd dips in 2001) and from CBI 
(2nd, 3rd, 4th, 5th peaks, 3rd, 4th dips at 1-sigma in 2002); (3) damping due to shear viscosity 
and the width of the region over which hydrogen recombination occurred when the universe was 
400000 years old (CBI 2002); (4) the primary anisotropies should have a Gaussian distribution (be 
maximally random) in almost all inflationary models, the best data on this coming from Boomerang; 
(5) secondary anisotropies associated with nonlinear phenomena subsequent to 400000 years, which 
must be there and may have been detected by CBI and another experiment, BIMA. Showing the 5 
"pillars" involves detailed confrontation of the experimental data with theory; e.g., (5) compares 
the CBI data with predictions from two of the largest cosmological hydrodynamics simulations 
ever done. DASI, Boomerang and CBI in 2002, AMiBA in 2003, and many other experiments 
have the sensitivity to demonstrate the next pillar, (6) polarization, which must be there at the 
~ 7% level. A broad-band DASI detection consistent with inflation models was just reported. A 
7th pillar, anisotropies induced by gravity wave quantum noise, could be too small to detect. A 
minimal inflation parameter set, {(O b , (Q cdm ,Q. tot ,Q.Q,WQ,n s , T c , <T 8 }, is used to illustrate the power 
of the current data. After marginalizing over the other cosmic and experimental variables, we find the 
current CMB+LSS+SN1 data give Q. tot = l.OO+'SL consistent with (non-baroque) inflation theory. 

Restricting to Q. tot = 1, we find a nearly scale invariant spectrum, n s — 0.97+ 9?. The CDM density, 

^cdm = ^«/mh 2 = -12- 01' an( ^ baryon density, co b = £2 & h 2 = .022+ ^, are in the expected range. 
(The Big Bang nucleosynthesis estimate is 0.019 ± 0.002.) Substantial dark (unclustered) energy 
is inferred, Q.q w 0.68 ± 0.05, and CMB+LSS i2g values are compatible with the independent 
SN1 estimates. The dark energy equation of state, crudely parameterized by a quintessence-field 
pressure-to-density ratio Wq, is not well determined by CMB+LSS (Wq < —0.4 at 95% CL), 
but when combined with SN1 the resulting Wq < —0.7 limit is quite consistent with the w-g=— 1 
cosmological constant case. 



A SYNOPSIS OF CMB EXPERIMENTS 



We are in the midst of a remarkable outpouring of results from the CMB that has seen 
major announcements in each of the last three years, with no sign of abatement in 
the pace as more experiments are scheduled to release analyses of their results. This 
paper is an update of [[!]] to take into account how the new data have improved the case 
for primordial acceleration, and for acceleration occurring now. The simplest inflation 
models are strongly preferred by the data. This does not mean inflation is proved, it 
just fits the available information better than ever. It also does not mean that competitor 
theories are ruled out, but they would have to look awfully like inflation for them to 
work. As a result many competitors have now fallen into extreme disfavour as the data 
have improved. 

The CMB Spectrum: The CMB is a nearly perfect blackbody of 2.725 ± 0.002^ [g], 
with a 3.372 ± 0.007 mK dipole associated with the 300 km s _1 flow of the earth in the 
CMB, and a rich pattern of higher multipole anisotropics at tens of /iK arising from 
fluctuations at photon decoupling and later. Spectral distortions from the blackbody 
have been detected in the COBE FIRAS and DIRBE data. These are associated with 
starbursting galaxies due to stellar and accretion disk radiation downshifted into the 
infrared by dust then redshifted into the submillimetre; they have energy about twice 
all that in optical light, about a tenth of a percent of that in the CMB. The spectrally 
well-defined Sunyaev-Zeldovich (SZ) distortion associated with Compton-upscattering 
of CMB photons from hot gas has not been observed in the spectrum. The FIRAS 95% 
CL upper limit of 6.0 x 10 5 of the energy in the CMB is compatible with the < 10 
expected from clusters, groups and filaments in structure formation models, and places 
strong constraints on the allowed amount of earlier energy injection, e.g., ruling out 
mostly hydrodynamic models of LSS. The SZ effect has been well observed at high 
resolution with very high signal-to-noise along lines-of-sight through dozens of clusters. 
The SZ effect in random fields may have been observed with the CBI and BIMA, again 
at high resolution, although multifrequency observations to differentiate the signal from 
the CMB primary and radio source contributions will be needed to show this. 

The Era of Upper Limits: The story of the experimental quest for spatial anisotropics 
in the CMB temperature is a heroic one. The original 1965 Penzias and Wilson discovery 
paper quoted angular anisotropics below 10%, but by the late sixties 10~ 3 limits were 
reached, by Partridge and Wilkinson and by Conklin and Bracewell. As calculations of 
baryon-dominated adiabatic and isocurvature models improved in the 70s and early 80s, 
the theoretical expectation was that the experimentalists just had to get to 10~ 4 , as they 
did, e.g., Boynton and Partridge in 73. The only signal found was the dipole, hinted at by 
Conklin and Bracewell in 73, but found definitively in Berkeley and Princeton balloon 
experiments in the late 70s, along with upper limits on the quadrupole. Throughout the 
1980s, the upper limits kept coming down, punctuated by a few experiments widely used 
by theorists to constrain models: the small angle 84 Uson and Wilkinson and 87 OVRO 
limits, the large angle 81 Melchiorri limit, early (87) limits from the large angle Tenerife 
experiment, the small angle RATAN-600 limits, the 7°-beam Relict- 1 satellite limit of 
87, and Lubin and Meinhold's 89 half-degree South Pole limit, marking a first assault 
on the peak. 

Primordial fluctuations from which structure would have grown in the Universe can be 



one of two modes: adiabatic scalar perturbations associated with gravitational curvature 
variations or isocurvature scalar perturbations, with no initial curvature variation, but 
variations in the relative amounts of matter of different types, e.g., in the number of 
photons per baryon, or per dark matter particle, i.e., in the entropy per particle. Both 
modes could be present at once, and in addition there may also be tensor perturbations 
associated with primordial gravitational radiation which can leave an imprint on the 
CMB. The statistical distribution of the primordial fluctuations determines the statistics 
of the radiation pattern; e.g., the distribution could be generically non-Gaussian, needing 
an infinity of N-point correlation functions to characterize it. The special case when 
only the 3D 2-point correlation function is needed is that of a Gaussian random field; 
if dependent upon spatial separations only, not absolute positions or orientations, it is 
homogeneous and isotropic. The radiation pattern is then a 2D Gaussian field fully 
specified by its correlation function, whose spherical transform is the angular power 
spectrum, where t is the multipole number. If the 3D correlation does not depend 
upon multiplication by a scale factor, it is scale invariant. This does not translate into 
scale invariance in the 2D radiation correlation, whose features reflect the physical 
transport processes of the radiation through photon decoupling. 

The upper limit experiments were in fact highly useful in ruling out broad ranges 
of theoretical possibilities. In particular adiabatic baryon-dominated models were ruled 
out. In the early 80s, universes dominated by dark matter relics of the hot Big Bang low- 
ered theoretical predictions by about an order of magnitude over those of the baryon- 
only models. In the 82 to mid-90s period, many groups developed codes to solve the 
perturbed Boltzmann-Einstein equations when such collisionless relic dark matter was 
present. Armed with these pre-COBE computations, plus the LSS information of the 
time, a number of otherwise interesting models fell victim to the data: scale invariant 
isocurvature cold dark matter models in 86, large regions of parameter space for isocur- 
vature baryon models in 87, inflation models with radically broken scale invariance lead- 
ing to enhanced power on large scales in 87-89, CDM models with a decaying (~ keV) 
neutrino if its lifetime was too long (> lOyr) in 87 and 91. Also in this period there were 
some limited constraints on "standard" CDM models, restricting €l tot , £l B , and the am- 
plitude parameter Og. (<7g is a bandpower for density fluctuations on a scale associated 
with rare clusters of galaxies, 8h~ ! Mpc, where h = // /(100 km s _1 Mpc -1 ).) 

DMR and Post-DMR Experiments to April 1999: The familiar motley pattern of 
anisotropics associated with 2 < i < 20 multipoles at the 30}lK level revealed by COBE 
at 7° resolution was shortly followed by detections, and a few upper limits (UL), at 
higher t in 19 other ground-based (gb) or balloon-borne (bb) experiments — most with 
many fewer resolution elements than the 600 or so for COBE. Some predated in design 
and even in data delivery the 1992 COBE announcement. We have the intermediate 
angle SP91 (gb), the large angle FIRS (bb), both with strong hints of detection before 
COBE, then, post-COBE, more Tenerife (gb), MAX (bb), MS AM (bb), white-dish (gb, 
UL), argo (bb), SP94 (gb), SK93-95 (gb), Python (gb), BAM (bb), CAT (gb), OVRO-22 
(gb), SuZIE (gb, UL), QMAP (bb), VIPER (gb) and Python V (gb). A list valid to April 
1999 with associated bandpowers is given in 0, and are referred here as 4.99 data. 
They showed evidence for a first peak although it was not well localized. A strong 
first peak, followed by a sequence of smaller peaks diminished by damping in the ^ 



spectrum was a long-standing prediction of adiabatic models. For restricted parameter 
sets, good constraints were given on n s , and on Q, tot and £l A when LSS was added [Q]. 

TOCO, BOOMERANG & MAXIMA: The picture dramatically improved over 
the 3 years since April 1999. In summer 99, the ground-based TOCO experiment in 
Chile [Q], and in November 99 the North American balloon test flight of Boomerang 
[P, gave results that greatly improved first-peak localization, pointing to Q. tot ~ 1. 
Then in April 2000 dramatic results from the first CMB long duration balloon (LDB) 
flight, Boomerang [0, |8|], were announced, followed in May 2000 by results from the 
night flight of Maxima [JTTJ|] . Boomerang's best resolution was 10.7' ± 1.4', about 40 
times better than that of COBE, with tens of thousands of resolution elements. (The 
corresponding Gaussian beam filtering scale in multipole space is £ s ~ 800.) Maxima 
had a similar resolution but covered an order of magnitude less sky. In April 2001, the 
Boomerang analysis was improved and much more of the data were included, delivering 
information on the spectrum up to £ ~ 1000 [ [TT| , |T2| ]. Maxima also increased its £ 
range [JT3p . 

Boomerang carried a 1.2m telescope with 16 bolometers cooled to 300 mK in the 
focal plane aloft from McMurdo Bay in Antarctica in late December 1998, circled the 
Pole for 10.6 days and landed just 50 km from the launch site, only slightly damaged. 
Maps at 90, 150 and 220 GHz showed the same basic spatial features and the intensities 
were shown to fall precisely on the CMB blackbody curve. The fourth frequency channel 
at 400 GHz is dust-dominated. Fig. [1] shows a 150 GHz map derived using four of the six 
bolometers at 150 GHz. There were 10 bolometers at the other frequencies. Although 
Boomerang altogether probed 1800 square degrees, the April 2000 analysis used only 
one channel and 440 sq. deg., and the April 2001 analyses used 4 channels and the 
region in the ellipse covering 800 sq. deg. That is the Boomerang data used in this paper. 
In [|T4"|], the coverage is extended to 1200 sq deg, 2.9% of the sky. 

Maxima covered a 124 square degree region of sky in the Northern Hemisphere. 
Though Maxima was not an LDB, it did well because its bolometers were cooled even 
more than Boomerang's, to 100 mK, leading to higher sensitivity per unit observing 
time, it had a star camera so the pointing was well determined, and, further, all frequency 
channels were used in its analysis. 

DASI, CBI & VSA: DASI (the Degree Angular Scale Interferometer), located at the 
South Pole, has 13 dishes of size 0.2m. Instead of bolometers, it uses HEMTs, operating 
at 10 frequency channels spanning the band 26 — 36GHz. An interferometer baseline 
directly translates into a Fourier mode on the sky. The dish spacing and operating 
frequency dictate the £ range. In DASI's case, the range covered is 125 < £ < 900. 
32 independent maps were constructed, each of size 3.4°, the field-of-view (fov). The 
total area covered was 288 sq. deg. DASI's spectacular results were also announced in 
April 2001, unveiling a spectrum close to that reported by Boomerang at the same time. 
The two results together reinforced each other and lent considerable confidence to the 
emerging ^ spectrum in the £ < 1000 regime. 

CBI (the Cosmic Background Imager), based at 16000 feet on a high plateau in 
Chile, is the sister experiment to DASI. It has 13 0.9m dishes operating in the same 
HEMT channels as DASI. The instrument measures 78 baselines simultaneously. The 
larger dishes by a factor of 4 and longer baselines imply higher resolution by about the 
same factor: the CBI results reported in May 2002 go to £ of 3500, a huge increase 
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FIGURE 1. The Boomerang 150 GHz bolometer map is shown in the top figure. Of the entire 1800 
square degrees covered, only the interior 800 sq. degs. (within the ellipse) were used in [|llj| and this 
analysis. In the April 2000 analysis only one channel within the rectangle was used. In the Ruhl et al., 
2002 analysis, essentially the entire region is covered. The three circles show regions cut out of the analysis 
because they contain quasars with emission at 150 GHz. The resolution is 10.7' ± lA'fwhm. 



over Boomerang, Maxima and DASI. Only the analyses of data from the year 2000 
observing campaign were reported. During 2000, CBI covered three deep fields of 
diameter roughly 0.75° [p)|], and three mosaic regions, each of size roughly 13 square 
degrees [JT7|] . In analyzing such high resolution data at 30 GHz, great attention must be 
paid to the contamination by point sources, but we are confident that this is handled well 
[[Rjp. Data from 2001 roughly doubles the amount, increases the area covered, and its 
analysis is currently underway. 

Fig. shows one of the CBI mosaic regions of the sky and Fig. [3| one of the deep 
regions. 

The VSA (Very Small Array) in Tenerife, also an interferometer, operating at 30 
GHz, covered the t range of DASI, and confirmed the spectrum emerging from the 
Boomerang, Maxima and DASI data in that region. The VSA is now observing at longer 
baselines to increase its i range. 

The Optimal Spectrum, circa Summer 2002: The power spectrum shown in Fig. || 
combines all of the data in a way that takes all of the uncertainties in each experi- 
ment (calibration and beam) into account. The point at small t is dominated by DMR, 
at 900 < I < 2000 by the CBI mosaic data, with that beyond 2000 by the CBI deep 
data. In between, Boomerang drives the small error bars, DASI and CBI set the cali- 
bration and beam of Boomerang and give spectra totally compatible with Boomerang. 
Both VSA and Maxima are in agreement with this data as well. Further, although 
the errors from the experiments before April 2000 are larger, the quite heterogeneous 
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FIGURE 2. The inner 2.5 degrees of one of the three 13 square degree CBI mosaic fields is shown. The 
mosaic image is a standard radio astronomy map, with point sources removed. The regions are smaller 
than Boomerang covered, but the resolution is a factor of at least three better. The mass subtended by the 
CBI resolution scale (~ 4') easily encompasses the mass that collapses later in the universe to generate 
clusters of galaxies. 



4.99+TOCO+Boomerang-NA mix of CMB data is very consistent with what the newer 
experiments show. It is an amazing concordance of data. Accompanying this story is a 
convergence with decreasing errors over time on the values of the cosmological param- 
eters given in Table |T]. 

Primary CMB Processes and Soundwave Maps at Decoupling: Boomerang, Max- 
ima, DASI, CBI and VSA were designed to measure the primary anisotropics of the 
CMB, those which can be calculated using linear perturbation theory. What we see in 
Figs. [1|, §are, basically, images of soundwave patterns that existed about 400,000 years 
after the Big Bang, when the photons were freed from the plasma. The visually evident 
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FIGURE 3. An 0.75 degree field-of-view (within the circle) Wiener-filtered map of one of the three 
deep CBI fields is shown. Apart from the primary anisotropy signal, the deep images contain power at 
high I that might be from the SZ effect in collapsed clusters at redshift of order 1. We know that the 
predicted SZ signal is near to that seen, and is very sensitive to a 8 C&j 52 ' « a§ ). <7 g w 1 is needed for the 
SZ effect as determined by hydro simulations to agree with the data. The primary CMB data prefer values 
between 0.8 and 0.9, and so although we do expect the SZ signal to be lurking within the CBI signal, it 
may not be quite as large as the extra power seen in these deep maps. 



structure on degree scales is even more apparent in the power spectra of the maps, which 
show a dominant (first acoustic) peak, with less prominent subsequent ones detected at 
varying levels of statistical significance. 

The images are actually a projected mixture of dominant and subdominant physical 
processes through the photon decoupling "surface", a fuzzy wall at redshift z d ~ 1050, 
when the Universe passed from optically thick to thin to Thomson scattering over a co- 
moving distance « 19 Mpc. Prior to this, acoustic wave patterns in the tightly-coupled 
photon-baryon fluid on scales below the comoving "sound crossing distance" at de- 
coupling, <, 150 Mpc (i.e., <, 150 kpc physical), were viscously damped, strongly so 
on scales below the ~ 10 Mpc damping scale. After, photons freely-streamed along 
geodesies to us, mapping (through the angular diameter distance relation) the post- 
decoupling spatial structures in the temperature to the angular patterns we observe now 
as the primary anisotropics. The maps are images projected through the fuzzy decou- 
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FIGURE 4. The ^ are defined in terms of CMB temperature anisotropy multipoles by ^ = £(£ + 
l)(|(Ar/r) ft?! | 2 )/(27r). The optimal ^ spectrum corresponds to a maximum-likelihood fit to the power in 
bands marginalized over beam and calibration uncertainties of the various experiments, This one uses "all- 
data" (DMR + Boomerang + Maxima + DASI + CBI mosaic + CBI deep + VSA + TOCO + Boomerang- 
NA + the 4.99 data). A A£ = 75 binning was chosen up to ~ 800, going over to the CBI deep binning at 
large i. The I > 2000 excess found with the deep CBI data is denoted by the light blue hatched region (95% 
confidence limit). Two best fit models to "all-data" are shown. They are both ACDM models; e.g., the 



upper curve (magenta) has parameters {£2 tot , £l t 



t c }={1.0,0.7, 0.02, 0.14,0.975,0}. 



pling surface of the acoustic waves (photon bunching), the electron flow (Doppler effect) 
and the gravitational potential peaks and troughs ("naive" Sachs-Wolfe effect) back then. 
Free-streaming along our (linearly perturbed) past light cone leaves the pattern largely 
unaffected, except that temporal evolution in the gravitational potential wells as the pho- 
tons propagate through them leaves a further AT imprint, called the integrated Sachs- 
Wolfe effect. Intense theoretical work over three decades has put accurate calculations 
of this linear cosmological radiative transfer on a firm footing, and there are speedy, 
publicly available and widely used codes for evaluation of anisotropics in a variety of 



cosmological scenarios, "CMBfast" and "CAMB" Q18|]. Extensions to more cosmologi- 
cal models have been added by a variety of researchers. 

Of course there are a number of nonlinear effects that are also present in the maps. 
These secondary anisotropics include weak-lensing by intervening mass, Thomson- 



scattering by the nonlinear flowing gas once it became "reionized" at z ~ 10, the thermal 
and kinematic SZ effects, and the red-shifted emission from dusty galaxies. They all 
leave non-Gaussian imprints on the CMB sky. 

The Immediate Future: Results are in and the analysis is underway for the bolometer 
single dish ACBAR experiment based at the South Pole (with about 5' resolution, 
allowing coverage to £ ~ 2000). Results from two balloon flights, Archeops and Tophat, 
are also expected soon. MINT, a Princeton interferometer, also has results soon to be 
released. 

In June 2001, NASA launched the all-sky HEMT-based MAP satellite, with 12' 
resolution. We expect spectacularly accurate results covering the £ range to about 600 to 
be announced in Jan. 2003, with higher £ and smaller errors expected as the observing 
period increases. Eventually four years of data are expected. 

Further downstream, in 2007, ESA will launch the bolometer+HEMT-based Planck 
satellite, with 5' resolution. 

Targeting Polarization of the Primary CMB Signal: The polarization dependence 
of Compton scattering induces a well defined polarization signal emerging from photon 
decoupling. Given the total ^ of Fig. |], we can forecast what the strength of that signal 
and its cross-correlation with the total anisotropy will be, and which £ range gives the 
maximum signal: ~ 5/lK over £ ~ 400 — 1600 is a target for the "E-mode" that scalar 
fluctuations give. 

A great race was on to first detect the it-mode. Experiments range from many degrees 
to subarcminute scales. DASI has analyzed 271 days of polarization data on 2 deep fields 
(3.4° fov) and just announced a 5a detection at a level consistent with inflation-based 
models [|T9[]. Boomerang will fly again, in December 2002, with polarization- sensitive 
bolometers of the kind that will also be used on Planck; as well, MAXIMA will fly again 
as the polarization-targeting MAXIPOL; the detectors on CBI have been reconfigured 
to target polarization and the CBI is beginning to take data; MAP can also measure 
polarization with its HEMTs. Other experiments, operating or planned include AMiBA, 
COMPASS, CUPMAP, PIQUE and its sequel, POLAR, Polarbear, Polatron, QUEST, 
Sport/BaRSport, BICEP, among others. The amplitude of the DASI E detection is 
0.8 ±0.3 of the forecasted amplitude from the total anisotropy data; the cross-correlation 
of the polarization with the total anisotropy has also been detected, with amplitude 
0.9 ± 0.4 of the forecast. These detections used a broad-band shape covering the £ 
range ~ 250 — 750 derived from the theoretical forecasts. The forecasts indicate solid 
detections with enough well-determined bandpowers for use in cosmological parameter 
studies are soon likely. 

We cannot yet forecast the strength of the "5-mode" signal induced by gravity waves, 
since there is as yet no evidence for or against them in the data. However, the amplitude 
would be very small indeed even at £ ~ 100. Nonetheless there are experiments such as 
BICEP (Caltech) being planned to go after polarization in these low £ ranges. 

Targeting Secondary Anisotropies: SZ anisotropics have been probed by single 
dishes, the OVRO and BIMA mm arrays, and the Ryle interferometer. Detections of 
individual clusters are now routine. The power at £ > 2000 seen in the CBI deep data 
(Fig. f§) and in the BIMA data at £ ~ 6000 may be due to the SZ effect in 
ambient fields, e.g., $23j\. A number of planned HEMT-based interferometers are being 
built with this ambient effect as a target: CARMA (OVRO+BIMA together), the SZA 



(Chicago, to be incorporated in CARMA), AMI (Britain, including the Ryle telescope), 
and AMiBA (Taiwan). Bolometer-based experiments will also be used to probe the SZ 
effect, including the CSO (Caltech submm observatory, a 10m dish) with BOLOCAM 
on Mauna Kea, ACBAR at the South Pole and the LMT (large mm telescope, with 
a 50 metre dish) in Mexico. As well, APEX, a 12 m German single dish, Kobyama, 
a 10 m Japanese single dish, and the 100m Green Bank telescope can all be brought 
to bear down on the SZ sky. Very large bolometer arrays with thousands of elements 
are planned: the South Pole Telescope (SPT, Chicago) and the Atacama Cosmology 
Telescope (ACT, Princeton), with resolution below 2', will be powerful probes of the SZ 
effect as well as of primary anisotropics. 

Anisotropics from dust emission from high redshift galaxies are being targeted by 
the JCMT with the SCUBA bolometer array, the OVRO mm interferometer, the CSO, 
the SMA (submm array) on Mauna Kea, the LMT, the ambitious US/ESO ALMA mm 
array in Chile, the LDB BLAST, and ESAs Herschel satellite. About 50% of the submm 
background has so far been identified with sources that SCUBA has found. 

The CMB Analysis Pipeline: Analyzing Boomerang and other single-dish experi- 
ments involves a pipeline, reviewed in []2~5|], that takes the timestream in each of the 
bolometer channels coming from the balloon plus information on where it is pointing 
and turns it into spatial maps for each frequency characterized by average temperature 
fluctuation values in each pixel (Fig. |T[) and a pixel-pixel correlation matrix characteriz- 
ing the noise. From this, various statistical quantities are derived, in particular the tem- 
perature power spectrum as a function of multipole, grouped into bands, and band-band 
error matrices which together determine the full likelihood distribution of the bandpow- 
ers p0| Bp. Fundamental to the first step is the extraction of the sky signal from the 
noise, using the only information we have, the pointing matrix mapping a bit in time 
onto a pixel position on the sky. In the April 2001 analysis of Boomerang, and subse- 
quent work, powerful use of Monte Carlo simulations was made to evaluate the power 
spectrum and other statistical indicators in maps with many more pixels than was possi- 
ble with conventional matrix methods for estimating power spectra JTT| ] . 

For interferometer experiments, the basic data are visibilities as a function of baseline 
and frequency, with contributions from random detector noise as well as from the sky 
signals. As mentioned above, a baseline is a direct probe of a given angular wavenumber 
vector on the sky, hence suggests we should make "generalized maps" in "momentum 
space" (i.e., Fourier transform space) rather than in position space, as for Boomerang. 
A major advance was made in [ PT| ] to deal with the large volume of interferometer data 
that we got with CBI, especially for the mosaics with their large number of overlapping 
fields. We "optimally" compressed the > <^(10 5 ) visibility measurements of each field 
into a < <^(10 4 ) coarse grained lattice in momentum space, and used the information in 
that "generalized pixel" basis to estimate the power spectrum and statistical distribution 
of the signals. 

There is generally another step in between the maps and the final power spectra, 
namely separating the multifrequency spatial maps into the physical components on 
the sky: the primary CMB, the thermal and kinematic Sunyaev-Zeldovich effects, the 
dust, synchrotron and bremsstrahlung Galactic signals, the extragalactic radio and sub- 
millimetre sources. The strong agreement among the Boomerang maps indicates that to 
first order we can ignore this step, but it has to be taken into account as the precision 



increases. 

Because of the 1 cm observing wavelength of CBI and its resolution, the contribution 
from extragalactic radio sources is significant. We project out of the data sets known 
point sources when estimating the primary anisotropy spectrum by using a number of 
constraint matrices. The positions are obtained from the (1 .4 GHz) NVSS catalog. When 
projecting out the sources we use large amplitudes which effectively marginalize over 
all affected modes. This insures robustness with respect to errors in the assumed fluxes 
of the sources. The residual contribution of sources below our known-source cutoff is 
treated as a white noise background with an amplitude (and error) estimated as well from 
the NVSS database flnj. 

The CMB Statistical Distributions are Nearly Gaussian: The primary CMB fluc- 
tuations are quite Gaussian, according to COBE, Maxima, and now Boomerang and 
CBI analyses. Analysis of data like that in the Fig. [I] map show a one-point distribution 
of temperature anisotropy values that is well fit by a Gaussian distribution. Higher or- 
der (concentration) statistics (3,4-point functions, etc.) tell us of non-Gaussian aspects, 
necessarily expected from the Galactic foreground and extragalactic source signals, but 
possible even in the early Universe fluctuations. For example, though non-Gaussianity 
occurs only in the more baroque inflation models of quantum noise, it is a necessary out- 
come of defect-driven models of structure formation. (Peaks compatible with Fig. |] do 
not appear in non-baroque defect models, which now appear highly unlikely.) There is 
currently no evidence for a breakdown of the Gaussianity in the 150 GHz maps as long 
as one does not include regions near the Galactic plane in the analysis. We have also 
found that the one-point distribution of the CBI data is also compatible with a Gaus- 
sian. However, since we know non-Gaussianity is necessarily there at some level, more 
exploration is needed. 

Though great strides have been made in the analysis of the Boomerang- style and CBI- 
style experiments, there is intense effort worldwide developing fast and accurate algo- 
rithms to deal with the looming megapixel datasets of LDBs and the satellites. Dealing 
more effectively with the various component signals and the statistical distribution of 
the errors resulting from the component separation is a high priority. 



COSMIC PARAMETER ESTIMATION 

Parameters of Structure Formation: Following [jlj], we adopt a restricted set 
of 8 cosmological parameters, augmenting the basic 7 used in [j^, 0, [FT], ^2p, 
{Q. A ,Q. k , 0) b ,0) cdm ,n s ,X c ,Og}, by one. The vacuum or dark energy encoded in the 
cosmological constant £l A is reinterpreted as Q.q, the energy in a scalar field Q which 
dominates at late times, which, unlike A, could have complex dynamics associated with 
it. Q is now often termed a quintessence field. One popular phenomenology is to add 
one more parameter, Wq = Pq/ Pq, where Pq and pQ are the pressure and density of the 

£>-field, related to its kinetic and potential energy by p Q = Q 2 /2 + (V£>) 2 /2 + V(Q), 

Pq = Q 2 /2 — (V<2) 2 /6 — V(Q). Thus w Q = — I for the cosmological constant. Spatial 
fluctuations of Q are expected to leave a direct imprint on the CMB for small i. This 
will depend in detail upon the specific model for Q. We ignore this complication here, 



but caution that using DMR data which is sensitive to low I behaviour in conjunction 
with the rest will give somewhat misleading results. To be self consistent, a model must 
be complete: e.g., even the ludicrous models with constant Wq would have necessary 
fluctuations to take into account. As well, as long as Wq is not exactly —1, it will vary 
with time, but the data will have to improve for there to be sensitivity to this, and for 
now we can just interpret Wq as an appropriate time-average of the equation of state. 
The curvature energy Q. k = 1 — Q. tot also can dominate at late times, as well as affecting 
the geometry. 

We use only 2 parameters to characterize the early universe primordial power spec- 
trum of gravitational potential fluctuations 4>, one giving the overall power spectrum am- 
plitude ^^(k n ), and one defining the shape, a spectral tilt n s (k n ) = 1 + dln&^/dlnk, 
at some (comoving) normalization wavenumber k n . We really need another 2, & GW {k n ) 
and n t {k n ), associated with the gravitational wave component. In inflation, the ampli- 
tude ratio is related to n t to lowest order, with &{n s — n t ) corrections at higher order, 
e.g., []2^|]. There are also useful limiting cases for the n s — n t relation. However, as one 
allows the baroqueness of the inflation models to increase, one can entertain essentially 
any power spectrum (fully ^-dependent n s {k) and n t {k)) if one is artful enough in de- 
signing inflaton potential surfaces. Actually n t (k) does not have as much freedom as 
n s (k) in inflation. For example, it is very difficult to get n t (k) to be positive. As well, one 
can have more types of modes present, e.g., scalar isocurvature modes (&»),«,■ Jk)) 
in addition to, or in place of, the scalar curvature modes (^ > ^(k n ),n s (k)). However, our 
philosophy is consider minimal models first, then see how progressive relaxation of the 
constraints on the inflation models, at the expense of increasing baroqueness, causes the 
parameter errors to open up. For example, with COBE-DMR and Boomerang, we can 
probe the GW contribution, but the data are not powerful enough to determine much. 
Planck can in principle probe the gravity wave contribution reasonably well. 

We use another 2 parameters to characterize the transport of the radiation through the 
era of photon decoupling, which is sensitive to the physical density of the various species 
of particles present then, ft). = ^ ; h 2 . We really need 4: a> b for the baryons, 0) cdm for the 
cold dark matter, 0£> hdm for the hot dark matter (massive but light neutrinos), and C0 er for 
the relativistic particles present at that time (photons, very light neutrinos, and possibly 
weakly interacting products of late time particle decays). For simplicity, though, we 
restrict ourselves to the conventional 3 species of relativistic neutrinos plus photons, with 
0) er therefore fixed by the CMB temperature and the relationship between the neutrino 
and photon temperatures determined by the extra photon entropy accompanying e + e 
annihilation. Of particular importance for the pattern of the radiation is the (comoving) 
distance sound can have influenced by recombination (at redshift z dec = cT, — 1), 

r, = 6000/^ Mpc^ adeC {(Om + co er a^y l/2 (l + co b a/(Aco r /3)y l/z d^/a, (1) 

where (Dy = 2.46 x 10~ 5 is the photon density, co er = 1.68ft)y for 3 species of massless 
neutrinos and co m = (O hdm + co cdm + co b . 

The angular diameter distance relation maps spatial structure at photon decoupling 
perpendicular to the line-of- sight with transverse wavenumber k ± to angular structure, 
through £ = & dec k,. In terms of the comoving distance to photon decoupling (recombi- 



nation), % dec , and the curvature scale d k , & dec is given by 



&dec = {ty^XdJ^'XdecAMZaec/dJ}, where d k = 3000|co,| l ' 2 Mpc, 
and x d ec = 6000 Mpc / (co m + (0 a~ 6w Q + co k ay l/2 dy/a. (2) 

The 3 cases are for negative, zero and positive mean curvature. Thus the mapping 
depends upon (O k , (Oq and Wq as well as on (O m . The location of the acoustic peaks £ pk ■ is 
proportional to the ratio of & dec to r s , hence depends upon co b through the sound speed as 
well. Thus £ pk . defines a functional relationship among these parameters, a degeneracy 
[ETA that would be exact except for the integrated Sachs- Wolfe effect, associated with 



the change of <I> with time if Q.q or Cl k is nonzero. (If <I> vanishes, the energy of photons 
coming into potential wells is the same as that coming out, and there is no net impact of 
the rippled light cone upon the observed AT.) 

Our 7th parameter is an astrophysical one, the Compton "optical depth" x c from 
a reionization redshift z reh to the present. It lowers ^ by exp(— 2t c ) at ts in the 
Boomerang/CBI regime. For typical models of hierarchical structure formation, we 
expect % c < 0.2. It is partly degenerate with a 8 and cannot quite be determined at this 
precision by CMB data now. 

The LSS also depends upon our parameter set: an important combination is the 
wavenumber of the horizon when the energy density in relativistic particles equals the 
energy density in nonrelativistic particles: k^ q ~ 5r _1 h _1 Mpc, where T fa Q, m hQ,~^ 2 . 
Instead of ^^{k n ) for the amplitude parameter, we often use at £ = 10 for CMB 
only, and <t| when LSS is added. When LSS is considered in this paper, it refers to 
constraints on T+ (n s — l)/2 and lncr| that are obtained by comparison with the data on 
galaxy clustering, cluster abundances and from weak lensing [ |Z3"| , ^p. At the current 



time, the constraints from <7 8 from lensing and cluster abundances are stronger thoes 
from r, although, with the wealth of data emerging from the Sloan Digital Sky Survey 
and the 2dF redshift survey, shape should soon deliver more powerful information than 
overall amplitude. However, in the future, weak lensing will allow amplitude and shape 
to be simultaneously constrained without the uncertainties associated with the biasing 
of the galaxy distribution wrt the mass that the redshift surveys must deal with. 

When we allow for freedom in co er , the abundance of primordial helium, tilts of tilts 
{dn, . x (k n ) / dink, ...) for 3 types of perturbations, the parameter count would be 17, 
and many more if we open up full theoretical freedom in spectral shapes. However, as 
we shall see, as of now only 4 combinations can be determined with 10% accuracy 
with the CMB. Thus choosing 8 is adequate for the present; 7 of these are discretely 
sampled[p8||, with generous boundaries, though for drawing cosmological conclusions 
we adopt a weak prior probability on the Hubble parameter and age: we restrict h to lie 
in the 0.45 to 0.9 range, and the age to be above 10 Gyr. 

Peaks, Dips and £l tot , Q.q and Wq: For given co m and co b , we show the lines of 
constant £ pk ■ & dec /r s in the Q. k -Q.Q plane for Wq=—\ in Fig. |5], and in the Wq-Q.q 
plane for Cl tot =\ in Fig. [7|, using the formulas given above and in p7|]. 



Our current best estimate Q22| ] of the peak locations £ pk ■, using all current CMB data 



and the flat+wk-/*+LSS prior, are 222 ±3, 537 ±6, 823 ±45, 1138 ±45, 1437 ±59, 
obtained by forming exp < \nl pk - >, where the average and variance of \n£ pk - are 
determined by integrating over the probability-weighted database described above. The 
interleaving dips are at 411 ±5, 682 ±48, 1023±44, 1319±51, 1653 ±48. With just the 
data prior to April 2000, the first peak value was 224 ± 25, showing how it has localized. 
A quadratic fit sliding over the data can be used to estimate the peak and dip positions 
in a model-independent way. It gives numbers in good accord with those given here 
[ |T2"1 , [17t ], but of course with larger error bars, in some cases only one-sigma detections 
for the higher peaks and dips. 

The critical spatial scale determining the positions of the peaks is r s , found by av- 
eraging over the model space probabilities to be 145 ± 2 Mpc comoving, thus about 
140 kpc as the physical sound horizon at decoupling. Converting peaks in fc-space into 
peaks in £-space is obscured by projection effects over the finite width of decoupling 
and the influence of sources other than sound oscillations such as the Doppler term. The 
conversion into peak locations in ^ gives £ pk ■ ~ fjjn& dec /r s , where the numerically 
estimated /■ factor is rs 0.75 for the first peak, approaching unity for higher ones. Dip 
locations are determined by replacing j by j + 1/2. These numbers accord reasonably 
well with the ensemble-averaged £ pk / dip j given above. 

The strength of the overall decline due to shear viscosity and the finite width of 
the region over which hydrogen recombination occurs can also be estimated, R D = 
10 ± 3 Mpc, i.e., 10 kpc back then, corresponding to an angular damping scale t D = 
1358 ±22. 

The constant I k • lines look rather similar to the contours in Figs. showing that 
the &ci ec /r s degeneracy plays a large role in determining the contours. The figures also 
show how adding other cosmological information such as H Q estimation can break the 
degeneracies. The contours hug the £l k = line more closely than the allowed £ pk . 
band does for the maximum probability values of C0 m and C0 b , because of the shift in 
the allowed £ pk ■ band as co m and (O b vary in this plane. The (O b dependence in r s would 
lead to a degeneracy with other parameters in terms of peak/dip positions. However, 
relative peak/dip heights are extremely significant for parameter estimation as well, and 
this breaks the degeneracy. For example, increasing co b beyond the nucleosynthesis (and 
CMB) estimate leads to a diminished height for the second peak that is not in accord 
with the data. 

Marginalized Estimates of our Basic 8 Parameters: Table [j] shows there are strong 
detections with only the CMB data for Q. tot , co b and n s in the minimal inflation-based 8 
parameter set. The ranges quoted are Bayesian 50% values and the errors are 1-sigma, 
obtained after projecting (marginalizing) overall other parameters. With "all-data", CO cdm 
begins to localize, but more so when LSS information is added. Indeed, even with just 
the COBE-DMR+LSS data, (O cdm is already localized. That Q.q is not well determined is 
a manifestation of the Q. k -Q.Q near-degeneracy discussed above, which is broken when 
LSS is added because the CMB -normalized ct 8 is quite different for open cf. pure Q- 
models. Supernova at high redshift give complementary information to the CMB, but 
with CMB+LSS (and the inflation-based paradigm) we do not need it: the CMB+SN1 
and CMB+LSS numbers are quite compatible. In our space, the Hubble parameter, 
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FIGURE 5. This shows lines of constant £. l in the Q. k -Q.Q plane (assuming Wq=— 1, i.e., a cosmo- 
logical constant) for the {©,„, C0 b } shown, near their most probable values. The data give £ k l — 222 ± 3. 
The higher peaks and dips have similar curves, scaled about the probable values listed above. Lowering 
CO h increases the sound speed, decreasing £ • ., and varying co m also shifts it. The 0.64 < h < 0.82 (heavier 
shading, H Q ) and 13 < age < 15 (lighter shading, t Q ) ranges and decelerations q = 0,-1/3,-1/2 are also 
noted. The sweeping back of the £ pk ■ curves into the closed models as £2 A is lowered shows that even if 
Q. tot =\, the phase space results in a ID projection onto the H tot axis that would be skewed to £l tot > 1. 
This plot explains much of the structure in the probability contour maps derived from the data, Fig. ^. 




h = (JV ,(n .h 2 )) 1 ' 2 , and the age of the Universe, t , are derived functions of the Q .h 2 : 
representative values are given in the Table caption. CMB+LSS does not currently give 
a useful constraint on Wq, though Wq < —0.7 with SN1. The values do not change very 
much if rather than the weak prior on h, we use 0.72 ± 0.08, the estimate from the 
Hubble key project [fTH]. Indeed just the CMB data plus this restricted range for h and 
the restriction to Q, tot = 1 results in a strong detection of Q.q. Allowing for a neutrino 
mass [ |32"| ] changes the value of Q.q downward as the mass increases, but not so much as 
to make it unnecessary. 

The Future, Forecasts for Parameter Eigenmodes: We can also forecast dramat- 
ically improved precision with future LDBs, ground-based single dishes and interfer- 
ometers, MAP and Planck. Because there are correlations among the physical vari- 
ables we wish to determine, including a number of near-degeneracies beyond that for 
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FIGURE 6. 1,2,3-sigma likelihood contour shadings for "all-data" and the weak-H+age prior prob- 
ability in the Q, k — £2g plane. The first interior lines are the 1,2-sigma ones when the LSS constraint is 
added, the most interior are the contours when the Hubble key project constraint is applied. The supernova 
contour shadings fl29[ ] are also plotted. Note that the contours are near the £l k = line, but also follow a 
weighted average of the £ k l ~ 220 lines. This approximate degeneracy implies £2g is poorly constrained 
for CMB-only, but it is broken when LSS is added, giving a solid SN1 -independent £2g "detection". When 
the Hubble key project constraint on H Q is added, partial breaking of this degeneracy occurs as well, as is 
evident from Fig. and from this figure. 



Q. k -Q.Q p7[], it is useful to disentangle them, by making combinations which diago- 
nalize the error correlation matrix, "parameter eigenmodes" [ |2~6| , |27| ]. For this exercise, 
we will add co hdm and n t to our parameter mix, but set Wq=—1, making 9. (The ratio 
^G\v(kn)/ ^^(kn) is treated as fixed by n t , a reasonably accurate inflation theory result.) 
The forecast for Boomerang based on the 800 sq. deg. patch with four 150 GHz bolome- 
ters used is 4 out of 9 linear combinations should be determined to ±0. 1 accuracy. This 
is indeed what was obtained in the full analysis of CMB only for Boomerang+DMR. 
The situation improves for the satellite experiments: for MAP, with 2 years of data, we 
forecast 6/9 combos to ±0.1 accuracy, 3/9 to ±0.01 accuracy; for Planck, 7/9 to ±0.1 
accuracy, 5/9 to ±0.01 accuracy. While we can expect systematic errors to loom as the 
real arbiter of accuracy, the clear forecast is for a very rosy decade of high precision 
CMB cosmology that we are now fully into. 
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FIGURE 7. Lines of constant £ • l in the Wq-Q.q quintessence plane (with Q. tot =l) are shown for the 



most probable values of {©,„, C0 h }. Lines of constant deceleration parameter q = (Q,„ 
H Q and age t Q in the ranges indicated are also shown. 
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TABLE 1. Cosmological parameter values and their 1-sigma er- 
rors are shown, determined after marginalizing over the other 6 cos- 
mological and the various experimental parameters, for "all-data" 
and the weak prior (0.45 < h < 0.9, age > 10 Gyr). The LSS prior 
was also designed to be weak. The detections are clearly very stable 
if extra "prior" probabilities for LSS and SN1 are included. Simi- 
lar tables are given in p2[]. If Q wt is varied, but Wq = — 1, param- 
eters derived from our basic 8 come out to be: age=15.2± 1.3 Gyr, 
h = 0.55 ± 0.08, £l m = 0.46 ± .1 1, Q b = 0.070 ± .02. Restriction to 
a tot = 1 and w Q = -1 yields: age= 14.1 ±0.6 Gyr, h = 0.65 ± .05, 
Q. m = 0.34 ± .05, £l b — 0.05 ± .006; allowing Wq to vary yields quite 
similar results. 
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